Probabilistic analysis of agent-based opinion formation models

When agent-based models are developed to capture opinion formation in large-scale populations, the opinion update equations often need to embed several complex psychological traits. The resulting models are more realistic, but also challenging to assess analytically, and hence numerical analysis techniques have an increasing importance in their study. Here, we propose the Qualitative Outcome Likelihood (QOL) analysis, a novel probabilistic analysis technique aimed to unravel behavioural patterns and properties of agent-based opinion formation models, and to characterise possible outcomes when only limited information is available. The QOL analysis reveals which qualitative categories of opinion distributions a model can produce, brings to light their relation to model features such as initial conditions, agent parameters and underlying digraph, and allows us to compare the behaviour of different opinion formation models. We exemplify the proposed technique by applying it to four opinion formation models: the classical Friedkin-Johnsen model and Bounded Confidence model, as well as the recently proposed Backfire Effect and Biased Assimilation model and Classification-based model.


Methods
The Qualitative Outcome Likelihood (QOL) analysis technique relies on QOL Tables and QOL Figures.We call QOL Tables the collection of five tables containing the probabilities that the model yields a predicted opinion distribution qualitatively categorised as PC, Co, Po, Cl, Di, given approximate information about the initial opinions and agent parameters.Conversely, QOL Figures are a collection of ten figures that contain the probabilities of obtaining a certain qualitative outcome, when only approximate information about the initial opinions or the agent parameters is known.
In order to construct QOL Tables and QOL Figures, we first need to compute the Qualitative Outcome Probability, i.e., the probability that the predicted opinion distribution belongs to a specific qualitative category.

Qualitative outcome probability
We consider a population with N agents, where agent i is associated with its time-varying opinion x i ∈ [−1, 1] , i ∈ V = {1, 2, . . ., N} , with x i [k] denoting the opinion of agent i at time k ∈ N .The opinion value represents the level of agreement or disagreement of the agent with a given statement: opinion values of +1 , −1 , and 0 represent complete agreement, complete disagreement, and indifference, respectively, while intermediate values correspond to less extreme opinions.The opinion distribution x is the set of all the opinions (non-ordered) of all the agents in the population.Each agent i ∈ V is also characterised by some parameters p i ; p denotes the set collecting all the parameter values (non-ordered) of all agents.
We then address the following question Q: for a population of N agents whose opinions evolve according to a given opinion formation model, if the initial opinion distribution belongs to a subset O of all possible initial opinion distributions, the agent parameters belong to a subset P of all possible agent parameters, and the underlying digraph belongs to a subset N of all possible digraphs with N vertices, what is the probability that, after K time steps, the predicted opinion distribution can be categorised as perfect consensus (or consensus, or polarisation, or clustering, or dissensus)?
This is a natural question to ask when only limited or incomplete information on the social system is available, which is the case in reality.When the information is not complete, the predicted opinions cannot be computed precisely, and one can at most evaluate the probability that they belong to a qualitative category.In some cases, knowing the exact initial opinions, agent parameters, and underlying digraph may not be very relevant (for instance, if the answer to Q with perfect consensus is 99% , then the qualitative outcome of the model is almost surely perfect consensus).
To provide a quantitative answer to question Q, we can proceed as follows: 1. Create three algorithms that randomly sample from the sets O , P , and N uniformly.2. For a number of events N e , do the following: (a) Randomly uniformly select an element from O , P , and N (an initial opinion distribution x, a set of agent parameters p, and an underlying digraph associated with weight matrix W), (b) evolve the system with the selected initial opinion distribution, agent parameter set, and underlying digraph according to the opinion model for the designated time steps, (c) categorise the predicted opinion distribution according to the sorting algorithm (based on 38 ; for details, see the Supplementary Information): if the resulting category is the desired one, then the event is labelled as a success.
3. The previous sequence of N e events has the structure of a Bernoulli process, where the random variable is a success if the predicted opinion distribution belongs to the chosen category.Hence, the number of events N e and successes N s can be used to suitably approximate the probability that the predicted opinion distribution actually belongs to the desired category, i.e., to answer question Q.
Then, the probability P B (O, P, N , K) that -if the initial opinion distribution belongs to the set O , the agent parameter set belongs to the set P , and the underlying digraph belongs to the set N -the predicted opinion distribution after K time steps is categorised as B ∈ {PC, Co, Po, Cl, Di} can be approximated by the Wilson score interval 43 : where z is the z-value for 95% confidence level, while N e and N s are respectively the number of events and of suc- cesses of a Bernoulli process where the Bernoulli trial is the answer to the question: Does the predicted opinion distribution belong to category B?
The higher the number of events, the better the approximation: in fact, δ → 0 as N e → ∞ .This procedure can be adjusted to account for different forms of available information.For instance, if some correlation between initial opinion assignation and agent parameters is known, it can be added as a constraint to the process.The only requirement is that the random variables are independent and identically distributed, which means that every possible initial configuration should be equally likely to be chosen so as to be evolved and produce predicted opinions that are then categorised.
(1) www.nature.com/scientificreports/In Eq. (1), P B is the real probability that answers question Q, while P is a value that we can compute to approximate the real probability with a desired confidence level.We choose the z-value for the 95% confidence level so that the probability P B (O, P, N , K) of the Bernoulli process belongs to the interval [P − δ, P + δ] with a probability of 0.95 (namely, the inclusion P B (O, P, N , K) ∈ [P − δ, P + δ] in Equation (1) is true with prob- ability 95% ).Other confidence levels may be chosen by using the corresponding z-value.
For all the simulations in this paper, the number of events N e is at least 10000.Hence, for any possible number of successes N s ∈ {0, 1, . . ., N e } , the uncertainty δ is less than 0.01 (see Fig. 1 in the Supplementary Information).Hence, for all the results we show, |P B − P| < 0.01 with 95% probability, where P is as defined in Eq. (1).For simplicity, we will shortly say that we are 'computing the probability P B ' when actually we are 'computing the value P such that |P B − P| < 0.01 with 95% probability' .
Equation ( 1) is useful to answer questions like Q, where the sets O , P , and N are given.However, a systematic application of Eq. (1) over a family of sets of initial opinion distributions, a family of sets of agent parameter sets, and a family of sets of underlying digraphs can provide further insight into the model properties and behaviour patterns.Here, we just present one possible application of this systematic approach, aimed at investigating the relation between P B and the initial opinion distribution and agent parameter set; but the approach could be adapted to investigate the relation between the model outcome and properties of the digraph, or other properties of the model.

Qualitative outcome likelihood tables and figures
Before analysing models through the systematic application of Eq. (1), we need to introduce some notation and concepts.

Families of opinion distribution sets and sets of agent parameter sets
Every opinion distribution x ∈ [−1, 1] N can be related to two values, the average x and the average of the absolute values |x| .Then, an opinion distribution can be represented as a point in the Cartesian plane, whose abscissa is |x| and whose ordinate is x , resulting in the Agreement Plot 39 .
We denote by {O i } an ordered family of opinion distribution sets.The generic set O i = O i (x) i , (|x|) i , ǫ o of opinion distributions is associated with three numbers, (x) i , (|x|) i and ǫ o > 0 , so that a generic opinion distribu- tion ξ belongs to O i if Agents are typically characterised by time-invariant parameters that capture their psychological traits.For instance, the generic agent i ∈ V is characterised by susceptibility i ∈ [0, 1] in the FJ model, by the triple the CB model, by entrenchment ρ i > 0 in the BEBA model, and by a confidence radius r ∈ [0, 2] in the BC model.The set of all agent parameters of a whole population is called the agent parameter set and generically denoted p (for the FJ, CB, BEBA, and BC models we denote the agent parameter set as , ψ , ρ , and r, respectively).
Every agent parameter set can also be related to two values, chosen depending on the specific model, which are suitably chosen functions of p that we generically denote as f a (p) and f b (p) .We consider: for the FJ model, the average and the variance σ ( ) of the susceptibility; for the CB model, the average radicalism β and average stubbornness γ ; for the BEBA model, the average ρ and the variance σ (ρ) of the entrenchment; for the BC model, the average r and the variance σ (r) of the confidence radius.Then, an agent parameter set can be represented as a point in the Cartesian plane, whose abscissa is f a (p) and whose ordinate is f b (p) , resulting in the Parameter Plot.
We denote by {P j } an ordered family of sets of agent parameter sets.The generic set P j = P j (p a,j , p b,j , ǫ p ) of agent parameter sets is associated with three numbers, p a,j , p b,j and ǫ p > 0 , so that a generic agent parameter set π belongs to P j if where the functions f a and f b are as previously described.

Implementation
Having defined the families {O i } and {P j } and how their elements relate to the points (x, |x|) and (p a , p b ) , we can now describe the procedure to create QOL Tables and QOL Figures.In the following, this systematic approach is explained resorting to the FJ model as a case study.
We focus on assessing how P B changes depending on the global features of the initial opinion distribution x and the agent parameter set p (which corresponds to the agent susceptibility for the FJ model).The set of underlying digraphs N is fixed, and includes 1000 digraphs with Small-World topology (for more details, see the Supplementary Information), and the evolution horizon is fixed to K = 50 .Then, the probability becomes a function of O and P only: P B = P B O, P .
We assess the qualitative behaviour of the FJ model by considering a family of initial opinion distribution sets that contains 318 different sets, {O i } 318 i=1 , and a family of sets of agent parameter sets that contains 371 different sets, {P j } 371 j=1 .In particular: • Every opinion distribution ξ ∈ O i meets the condition in Eq. ( 2) with ǫ o = 0.05 .Figure 1a shows a repre- sentation of the 318 considered sets O i , where each green dot corresponds to a set O i , i = 1, . . ., 318 , and its ( location in the Cartesian plane is determined, based on Eq. ( 2), by the corresponding values of (x) i and (|x|) i .• Every set of agent parameters π ∈ P j with average π and variance σ (π) meets the condition in Eq. ( 3) with ǫ p = 0.05 , p a,j = ( ) j , f a (π) = π , p b,j = (σ ( )) j , and f b (π) = σ (π) .Figure 1b shows a representation of the 371 considered sets P j , where each magenta dot corresponds to a set P j , j = 1, . . ., 371 , and its location in the Cartesian plane is determined, based on Eq. ( 3), by the corresponding values of ( ) j and (σ ( )) j .
Thanks to this characterisation, we can denote with the probability that, starting from an initial opinion distribution with average x and average absolute value |x| , for a population whose susceptibility has average and variance σ ( ) and whose opinions evolve according to the FJ model over a strongly connected Small-World network digraph in N , the predicted opinion distribution after K = 50 time steps belongs to the category B , with B ∈ {PC, Co, Po, Cl, Di}.Any pair of points (one green and one magenta) in Fig. 1a and b represents a pair of sets O i and P j , which can be associated with a number in the interval [0, 1] that quantifies the probability P B = P B (O i , P j ) = P B (x) i , (|x|) i , ( ) j , (σ ( )) j .Since P B depends on four parameters, it cannot be plotted directly.
However, as shown by the schematic example in Fig. 2, by suitably ordering the points in Fig. 1a and b, all the probabilities P B for different choices of (x, |x|) and of ( , σ ( )) can be arranged in a Table .Since B ∈ {PC, Co, Po, Cl, Di} , this procedure actually leads to five different tables, one per qualitative category of predicted opinion distribution, which are collectively called the Qualitative Outcome Likelihood Tables (QOL Tables).
In the FJ model case, each of the five QOL Tables has 318 rows (the number of points in Fig. 1a) and 371 columns (the number of points in Fig. 1b), and each of its 117978 cells contains a number between 0 and 1. Due to the huge number of cells, instead of showing the value of each cell, it is more practical and insightful to represent the tables visually by an image, where the colour of each cell depends on its value: probabilities near 0 (respectively, 0.5, and 1) have a light blue (resp.purple, and orange) colour.This colour convention for the probability values is used for all the QOL analysis and can be seen in the colour-bars of the figures.
The order of the points in Fig. 1a and b is directly related to the order of rows and columns in the QOL Tables.Hence, a suitable order may reveal patterns, such as clusters or regions in the table with particularly high or low probability.One can also plot the histogram of the cell values in the QOL Tables, which is invariant with respect to row and column permutations and represents the overall probability that a certain predicted category is obtained.For the FJ model, Fig. 3a shows the order of points, while Fig. 4a shows the corresponding visualisations of the QOL Tables and their corresponding histograms.
Assume now that no information at all is available about the initial opinion distribution (respectively, about the agent parameters).Then, we can perform the analysis by taking O (resp.P ) as the set of all possible initial opinion distributions (resp.all possible agent parameter sets).In this case, P B depends only on and σ ( ) (resp. on x and |x| ), and therefore the resulting probabilities can be visualised directly in the Parameter Plot (resp.Agreement Plot).The resulting ten QOL Figures show the probabilities P B for B ∈ {PC, Co, Po, Cl, Di} when only incomplete information on either the initial opinions or the agent parameters is available.Figure 5a shows the results for the FJ model.QOL Figures can be seen as a visualisation of the QOL Tables when there is only one set of initial opinion distributions (respectively, of agent parameter sets): since there is only one row (resp.  with coordinates ( , σ ( )) , we can assess the probability P B that a population starting from an initial opinion distribution with average x and average of the absolute values |x| , whose susceptibility has average and variance σ ( ) , produces a predicted opinion distribution categorised as B , with B ∈ {PC, Co, Po, Cl, Di}.column), the table values can be directly represented in the Parameter Plot (resp.Agreement Plot).The resulting figures are invariant to the order of opinion distribution sets and of sets of agent parameter sets.
To summarise, the complete QOL analysis involves a three-step process: 1. Agreement Plot and Parameter Plot Ordering: To analyse the model, define the family of possible initial opinion distribution sets {O i } , where each set O i = O i (x) i , (|x|) i , ǫ o is characterised according to Eq. ( 2), and the family of possible sets of agent parameter sets {P j } , where each set P j = P j (p a,j , p b,j , ǫ p ) is charac- terised according to Eq. ( 3).Order them suitably (a suitable order is one that can be easily interpreted and that results in QOL Tables with clear regions of high probability; in this paper, the ordering was found by numerically exploring various possibilities).2. Compute and analyse the QOL Tables: For each pair of opinion distribution and agent parameter sets (O i , P j ) , compute the probability P B (O i , P j ) that the predicted opinion distribution is categorised as B for B ∈ {PC, Co, Po, Cl, Di} ; the probability in Eq. ( 4) is approximated by the value P in Eq. ( 1).This provides the value associated with row i and column j in the corresponding table.Then, plot the table by colour-coding the cells depending on their value and analyse them taking into account the row and column order interpretation.It is also possible to plot and analyse the table histograms, which are independent of the order chosen at the previous step 1.

Compute and analyse the QOL Figures: If no information on the initial opinion distribution is available,
the probability P B depends only on the considered set of agent parameter sets P j ; therefore, the probability values can be colour-coded in the Parameter Plot of the family {P j } .Analogously, if no information on the agent parameters is available, the probability P B depends only on the considered set of opinion distributions O i ; therefore, the probability values can be colour-coded in the Agreement Plot of the family {O i }.
If no information is available on both the initial opinions and the agent parameters, then one can only consider O as the set of all possible opinion distributions and P as the set of all possible agent parameter sets, which yields QOL Tables with a single cell, i.e. a single probability, for each table.Table 1 reports the five probabilities for each of the considered models.The procedure described above requires that the opinion x i of each agent is scalar-valued and belongs to the interval [−1, 1] .However, our analysis can be adapted to models for which this is not the case.If the opinions are still unidimensional and bounded (for instance, in the interval [0, 1] instead), our proposed procedure can be applied by either normalising the opinions to be in the [−1, 1] interval or suitably adapting the definition of the considered opinion categories (Perfect Consensus, Consensus, Polarisation, Clustering, and Dissensus).If the opinions are not bounded, or are multidimensional, then the procedure can be applied after having suitably ) corresponding to a set P j : they are ordered as P 1 , P 2 , P 3 , P 4 following the red arrows.Then, the QOL Tables are constructed, whose rows are associated with points in the Agreement Plot and whose columns are associated with points in the Parameter Plot, in the corresponding order.Therefore, the cell in position (i, j) in the QOL Tables is associated with point O i = (x) i , (|x|) i in the Agreement Plot and with point P j = (p a,j , p b,j ) in the Parameter Plot, and contains a value between 0 and 1 that represents the probability P B .
redefined the opinion categories.In models with multidimensional opinions, where all dimensions are bounded, it is also possible to perform our proposed probabilistic analysis in each dimension independently.
The procedure to compute QOL Tables and QOL Figures is summarised in the following pseudo-code: 1. Input: the number of agents N ; the simulation horizon K; a collection ϒ = x i , |x| i I i=1 representing I points in the Agreement Plane; a collection � = (p a ) j , (p b ) j J j=1 representing J points in the Parameter Plane; a collection of interaction networks of suitable size N ; a number of events N e ; tolerance values ǫ o and ǫ p ; and a colour map CM for the interval [0, 1]. 2. Initialise five empty tables T PC , T Co , T Po , T Cl , T Di with I rows and J columns. 3.For each point x i , |x| i in ϒ compute a set of initial opinion distributions O i satisfying Equation (2): the collection of these sets forms the ordered family {O i } I i=1 .4. For each point (p a ) j , (p b ) j in compute a set of agent parameter sets P j satisfying Equation (3): the col- lection of these sets forms the ordered family {P j } J j=1 .
5. For each pair of indices (i, j) where i = 1, . . ., I and j = 1, . . ., J do the following: (a) Initialise five counters N s,i,j,PC = 0 , N s,i,j,Co = 0 , N s,i,j,Po = 0 , N s,i,j,Cl = 0 , N s,i,j,Di = 0 : these are the success counters for each of the Bernoulli trials for the pair (i, j).(b) For N e iterations do: (i) Take a randomly selected opinion distribution ξ ∈ O i , a randomly selected set of agent parameters π ∈ P j , and a randomly selected underlying digraph G ∈ N .rows and low columns, which for the FJ model, according to the order in Fig. 3a, are associated with societies with high x and low |x| (i.e., initially polarised), and with high and low σ ( ) (i.e., very susceptible).
(ii) Evolve the model with agent parameters π , initial opinions ξ , over the digraph G for K time steps.Call the predicted opinions ξ f .(iii) Classify the predicted opinions ξ f into one of the five possible categories Perfect Consensus (PC), Consensus (Co), Polarisation (Po), Clustering (Cl), and Dissensus (Di).(iv) For B ∈ {PC, Co, Po, Cl, Di}: (A) Answer the Bernoulli trial: Does the predicted opinion distribution belong to category B? (B) If the answer is yes, increase the counter N s,i,j,B by one (iii) For B ∈ {PC, Co, Po, Cl, Di} (i) Compute the probability (ii) Store the probability P i,j,B in row i, column j of the table T B .
6. Plot the completed tables T PC , T Co , T Po , T Cl , T Di (each with I rows and J columns with values in the interval [0, 1]) using the colour map CM to obtain QOL Tables.If desired, also plot histograms of the values contained in the tables (the plots will be similar to the ones presented in Fig. 4).The rows corresponds to points (x, |x|) in the Agreement Plot and the columns to points (p a , p b ) in the Parameter Plot, both ordered as in Fig. 3: each cell is associated with a set of initial opinion distributions and a set of agent parameter sets.The cell colour represents the probability P B , computed as per Eq. 1 (with 10000 samples), that the predicted opinion distribution after 50 time steps is categorised as B ∈ {PC, Co, Po, Cl, Di} .For details on the construction of QOL Tables, see Fig. 2. Second row: histograms corresponding to the QOL Tables, visualising the number of cells in the table whose probability value P B is in a given interval; the bar colour reflects the corresponding probability.
Remark: if in the table plots there are no distinguishable patterns, it is possible to reorder the elements in the ordered families {O i } I i=1 and {P j } J j=1 (which is equivalent to permuting rows and columns of the table) so as to reveal hidden patterns.7. Since each element in {O i } I i=1 (respectively, {P j } J j=1 ) is associated with a point in the Agreement Plot (resp.Parameter Plot) by Eq. (2) (resp.(3)), the ordering can be represented by colour coding points in the Agreement Plot (resp.Parameter Plot).This will be useful in the analysis of the QOL Tables (resulting in plots such as the ones presented in Fig. 3). Figure 5.When no information on agent parameters (respectively, initial opinions) is available, the probabilities P B that the predicted opinion distribution after 50 time steps is categorised as B ∈ {PC, Co, Po, Cl, Di} can be visualised directly in the Agreement Plot (respectively, in the Parameter Plot), as shown in the first (respectively, second) row of each panel.For each dot, the probability P B is computed as per Eq. ( 1) with a different number of samples.The number of samples for each point in the first row is 10,000 times the number of Agent Parameter points shown in Fig. 3 (in the right plots of each panel).The number of samples for each point in the second row is 3,180,000, since that is 10,000 times the number of points in the Agreement Plot in Fig. 3 (in the left plots of each panel).(ii) Compute the probability (iii) Plot a marker in the Cartesian Plane in coordinates x i , |x| i .The colour of this marker is given by the value P i,•,B and the colour map CM.
(iv) For j = 1, . . ., J (i) Add all the Bernoulli trial successes of simulations for which the agent parameter set belonged to P j .Numerically, this is equivalent to: (ii) Compute the probability (iii) Plot a marker in the Cartesian Plane in coordinates (p a ) j , (p b ) j .The colour of this marker is given by the value P •,j,B and the colour map CM.
The results of the previous plots are equivalent to the ones presented in Fig. 5.
In principle, the Agreement Plot and Agent Parameters can be ordered as desired so as to produce QOL Tables where patterns are more clearly visible.However, if the ordering is too intricate, the interpretation of QOL Tables becomes difficult.This leads to a trade-off between producing easy-to-interpret orderings and QOL Tables with clear patterns.

Results
Following our proposed three-step process, described in detail in the Methods section, we analyse four opinion formation models: the Friedkin-Johnsen (FJ), Classification-based (CB), Backfire Effect and Biased Assimilation (BEBA), and Bounded Confidence (BC) models.

Agreement plot and parameter plot ordering
Figure 3 shows how the points in the Agreement Plot and in the Parameter Plot can be ordered (differently) for the FJ, CB, BEBA, and BC models.
For the FJ model, when the Agreement Plot points are ordered as in Fig. 3a (thus matching the order of rows in the QOL Tables), the numbering grows from the corners (0, 0) and (1, 1) to the corner (1, 0).Hence, the first rows of the corresponding QOL Tables are approximately associated with opinion distributions that are initially categorised as perfect consensus, in the following rows as consensus, then as clustering, then as dissensus, and eventually as polarisation in the last rows.As for the order of the Parameter Plot points, the first columns correspond to societies with very high susceptibility values, the following columns to societies with a lower susceptibility, while the last columns correspond to societies where every agent has the same value of susceptibility (around 0.3).The order is more strongly affected by the average than by the variance.
For the CB model, Fig. 3b shows that the first rows of the QOL Tables correspond to societies with very extreme opinions (high |x| ), while the last rows to societies that are almost completely indifferent.The first columns are associated with highly radical societies (high β ), while the last columns have on average a small radical weight and any proportion of conformist and stubborn weights.
For the BEBA model, the order of the Agreement Plot points is shown in Fig. 3c (and is quite similar to the one in Fig. 3a): the first rows in the QOL Tables approximately correspond to initial opinion distributions forming perfect consensus, while moving towards the last rows the initial opinion distribution changes from perfect consensus to consensus, dissensus, clustering and, finally, polarisation.On the other hand, the order of the Parameter Plot points is peculiar, and highly depends on both the average and the variance of the agents' entrenchment.The first columns correspond to low values of entrenchment and as the order increases two groups appear, one with high average entrenchment and the other with intermediate-high average entrenchment and high variance.Therefore, when analysing the QOL Tables for the BEBA model, we need to take into account that right columns may correspond to societies where the average entrenchment is very high, or where only some agents have a very high entrenchment (this would explain the average of around 4 and high variance).
For the BC Model, the Agreement Plot ordering of Fig. 3d means that the first rows in the QOL Tables correspond to perfect consensus and move through consensus, clustering, and dissensus until reaching polarisation for the last rows.The agent parameter ordering shows a significant dependence on the confidence radius variance: the first columns correspond to high mean radius with low variance and the last columns to intermediate-low mean confidence radius and high variance.Hence, in the societies represented in the left columns in the QOL Tables, most agents have the same confidence radius, around 2, and in societies to the right some agents have a high and others have a low confidence radius.

Qualitative outcome likelihood tables
Figure 4 presents the QOL Tables for the FJ, CB, BEBA, and BC models.Here we briefly discuss the results achieved for each of the models and highlight the peculiarities of their behaviour.

FJ model
Figure 4a shows that perfect consensus is an almost certain outcome for the first rows.In fact, since the opinion distributions in these rows are initially located around the corners (0, 0) and (1, 1), they start as perfect consensus and, given that the FJ model has no mechanism to diversify opinions, the opinion distributions remain categorised as perfect consensus.Perfect consensus is also almost sure for cells in the left columns, corresponding to societies with very high susceptibility.In fact, if all agents are highly susceptible, they reach perfect consensus, because in these cases the FJ model boils down to the French-DeGroot model.
Figure 4a also shows that polarisation can be expected only when agents start from an already polarised opinion distribution and the susceptibility is very low and uniform (see the orange bottom right corner in the third table): the agents remain with their initial opinion, which was polarised from the start.Also other opinion distribution categories have regions in which they are more likely, but with a less intuitive interpretation.
The histograms in Fig. 4a provide additional information that is independent of the row and column order.The only category which can be almost guaranteed as an outcome is perfect consensus: only in this case, the right bin (in the neighbourhood of probability 1) has a significant height.Also, the least likely outcome is polarisation: the corresponding histogram has the largest left bin (in the neighbourhood of probability 0).The consensus histogram has an almost uniform probability distribution, meaning that the outcome probability highly depends on missing information.

CB model
Only perfect consensus and consensus are likely to occur in well-defined regions in the QOL Tables in Fig. 4b, suggesting that this model's behaviour is highly intricate.Other Agreement Plot and Parameter Plot orders may create better defined regions, but would not follow an easy-to-interpret pattern and, as such, they would not provide interpretable information about the model behaviour.
The histograms show that, except for dissensus, all qualitative outcomes have a non-negligible probability of being achieved.For every category, a significant number of cells has a probability different from zero or one.Interestingly, perfect consensus and consensus are likely outcomes in opposite scenarios: perfect consensus occurs in the left columns, corresponding to highly radical societies, while consensus can be found in societies with low radicalism.This happens because radicalism moves opinions to extremes, which often results in perfect consensus (the only other outcome produced by radicalism is polarisation, but it requires the initial opinions to have almost null bias, which is rarely the case, hence the most likely outcome is perfect consensus, as one extreme will have significantly more agents that the other), while consensus requires a clear but not extreme tendency of all agents towards similar opinions, which is more likely achieved by high conformist traits.Clustering is partially present in the left side of the last rows, corresponding to significantly radical societies starting from almost indifferent opinions.In fact, individuals with radical traits may move from indifference towards a stronger opinion, but not to such an extent to produce perfect consensus: they may end up forming two or three not very distant subgroups that correspond to clustering.
Comparing Fig. 4a and b reveals that perfect consensus and consensus, which have a high probability in adjacent cells for the FJ model, have a high probability in opposite cells for the CB model.

BEBA model
Figure 4c shows that the BEBA model has a strong tendency to transform opinion distributions into either perfect consensus or polarisation.In particular, perfect consensus is by far the most likely outcome, which occurs almost always; in the rare cases in which perfect consensus is not achieved, then the outcome is polarisation.This can be somewhat unexpected, as backfire effect and biased assimilation are generally associated with polarisation.However, according to our definition, polarisation requires the existence of two 'distant enough' subgroups having 'comparable' size (The formal definition of 'distant enough' and 'comparable' group size depends on the interpretation and code implementation; the quantitative metrics adopted in the manuscript are explained in detail in the Supplementary Information file, where the algorithm to sort opinion distributions is described.In particular, in our case, 'distant enough' means that the difference between opinions is at least 0.6, while 'comparable size' means that the difference between the numbers of individuals in the two groups is small enough, at most amounting to 38% of the population in the most extreme case: in fact, none of the two groups includes more than 50% of the population or less than 12% of the population, but the two groups together include more than 50% of the population.).If this last condition is not met, e.g. because one group has significantly more individuals than the other, then the opinion distribution is categorised as perfect consensus.Although the BEBA model promotes www.nature.com/scientificreports/ the existence of these two separate subgroups, it includes no intrinsic mechanism to equalise the number of individuals in each one and, since random initial assignations most likely produce unequal subgroups, most of them are categorised as perfect consensus.
The conditions for which the BEBA model most likely produces polarisation can be inferred from Fig. 4c: polarisation mostly occurs in the lower right corner, which corresponds to initial opinions near the (1, 0) point in the Agreement Plot (associated with a polarised distribution), and agent parameters with either high average entrenchment or intermediate-high average entrenchment and high variance.Then, if the population is already polarised, high entrenchment ensures that the two separate subgroups forming the initial polarisation distribution are maintained.

BC model
The QOL Tables in Fig. 4d show that the Bounded Confidence model predominantly produces perfect consensus and that the second most likely outcome is clustering.Perfect consensus is almost guaranteed for the left part of the QOL Tables, that is, approximately when the mean confidence radius is above 1.This makes sense: with such a high mean confidence radius, even with high variance, several agents will be influenced by almost all the other agents and cause global convergence to a single opinion, producing perfect consensus.
The QOL Tables show that clustering is almost guaranteed for the bottom right corner, that is, when the initial opinions are near the point (1, 0) (most likely polarised), and the confidence radius has an intermediate-low mean (around 0.5) and high variance (around 0.8).This happens because, when the initial opinions start from a polarised condition, the agents with intermediate confidence radius move towards the agents with low confidence radius, which remain in place, as they are influenced by very few agents.This behaviour has the potential of moving agents away from the two initial polarised groups and thus lead to clustering.
From Fig. 4d it is also clear that the other qualitative outcomes are extremely unlikely, because, in order to exist, they would need to start from very specific initial opinions, and the initial opinions and agent parameters would require a very specific configuration, which is not likely to happen at random.

FJ model
Figure 5a confirms that perfect consensus is the only outcome that occurs with very high probability for some regions of the Agreement Plot and the Parameter Plot.Looking at the Agreement Plot, this happens when the average of the opinions and of their absolute values is in the the neighbourhood of corners (0, 0) or (1, 1): opinions starting from perfect consensus are very likely to remain in that category, since the FJ model includes no mechanism to move opinions away from each other.Also, looking at the Parameter Plot, perfect consensus is very likely when the average susceptibility is high: in this case, the model behaves as the French-DeGroot model, which is known to lead to perfect consensus.
In addition to these clear trends, analysing the other plots of Fig. 5a reveals interesting observations.For both the Agreement Plot and the Parameter Plot, the consensus probability is significant (around 50%) for a wide variety of initial opinions and agent parameters, as shown by the vast purple area.Conversely, the probabilities of obtaining other opinion categories are low.Polarisation can only occur in the FJ model when the opinions are initially polarised and the susceptibility is low.Similarly, clustering can only happen when the opinions start either in clusters or dissensus and the agent parameters are 'uniformly' allocated, so that in the proximity of different opinions some agents have a low susceptibility (and thus keep their own opinion) while others have high susceptibility (and thus move towards the agents with low susceptibility).Although this combination of circumstances is unlikely, it is still noticeable in the plots corresponding to clustering, especially in the purple colour of Parameter Plot points associated with high variance and low average of the susceptibility.

CB model
Unlike the QOL Tables in Fig. 4b, the plots in Fig. 5b show clear regions in the Agreement Plot and Parameter Plot where specific qualitative outcomes are likely to be achieved.Perfect consensus, polarisation and clustering have clear regions in the Agreement Plot associated with significant probability: the neighbourhood of points (1, 1), (1, 0), and (0, 0) for perfect consensus, polarisation, and clustering respectively.For perfect consensus and polarisation, those opinions already belong to the desired qualitative category, therefore, if they do not change much, their qualitative category will not change.The possible explanation for the clustering case was previously discussed.
In the Parameter Plot, perfect consensus, consensus, polarisation, and clustering have regions associated with significant probability, which seem to depend mostly on the average radical weight.From minimum to maximum radical weight, the order in which the qualitative outcomes are more likely is consensus, clustering, polarisation, and perfect consensus.In fact, the existence of these categories requires a degree of intermediate and not extreme opinions: consensus needs multiple similar but not identical opinions, clustering needs two or more subgroups with distinguishable opinions, polarisation needs two subgroups with a comparable number of agents having distant opinions, and perfect consensus needs most agents having very similar, and possibly extreme, opinions.The average radical weight directly relates to how extreme opinions are in the population.

BEBA model
Figure 5c (row 2, column 1) shows that perfect consensus is by far the most likely outcome, even in the presence of high entrenchment.This further corroborates the conclusion, drawn from Fig. 4c, that high entrenchment does not necessarily lead to polarisation.The Parameter Plots in Fig. 5c show that, when all initial opinions are considered, only intermediate-low average values of entrenchment with high variance have a significant probability of producing opinions that are not categorised as perfect consensus.

BC model
Figure 5d provides more specific information about the conditions for which perfect consensus is not the almost sure outcome.This happens when the variance is maximal and the mean confidence radius is below 1 (see the purple area in row 2, column 1 of Fig. 5d).In that region, there is also a significant probability of obtaining clustering, and when the mean confidence radius is below 0.5 there is a slight probability of obtaining polarisation.This happens, however, not because opinion distributions from different categories evolve towards polarisation, but because polarised opinions remain polarised.

Global model probability outcomes
Table 1 reports the outcome probabilities when no information at all is known about both the initial opinion distribution and the agent parameters, and therefore shows the approximate probabilities that an opinion formation model yields each of the five possible outcomes.Again the FJ, CB, BEBA, and BC models are considered.Interestingly, the FJ model is the only for which perfect consensus is not the most likely outcome, despite being one of the two models without a mechanism to disperse opinions (the other one is the BC model): in fact, the susceptibility trait in this model allows opinions to move close, but prevents them from converging to the exact same opinion.Both the FJ and the CB models show a remarkable ability to yield all the different qualitative behaviours, although polarisation is relatively unlikely for the FJ model and dissensus is relatively unlikely for the CB model.Conversely, the BC and BEBA models are strongly biased towards perfect consensus.It can also be observed that, for the BEBA and BC models, the second most likely outcome (polarisation for the BEBA model, clustering for the BC model) has approximately the same probability.

Discussion
We have proposed a novel numerical simulation-based analysis technique for agent-based opinion formation models.Starting from incomplete information on the considered population (initial opinion distribution, agent parameters, underlying digraph), we can approximate the probability that the predicted opinion distribution belongs to a qualitative category (perfect consensus, consensus, polarisation, clustering, or dissensus).The systematic application of this approximation can be used to construct Qualitative Outcome Likelihood Tables and Figures, which reveal intrinsic model properties and behaviour patterns.We have applied the proposed technique to study four agent-based opinion formation models: Friedkin-Johnsen (FJ), Classification-based (CB), Backfire Effect and Biased Assimilation (BEBA), and Bounded Confidence (BC).
For the FJ model, our results suggest that the average population susceptibility is a good indicator of the predicted opinion distribution category.Moving from complete average susceptibility = 1 to null average sus- ceptibility = 0 , the most likely qualitative outcomes are perfect consensus, consensus, dissensus, clustering, and polarisation.Polarisation is the least likely outcome and only occurs when the initial condition is already at polarisation: in fact, the FJ model includes no mechanism to move opinions apart.For the same reason, the only outcome that occurs almost certainly in some conditions (when the initial opinions are either indifferent, near 0, or extreme, near 1 or −1 ) is perfect consensus.
However, including a mechanism that drives opinions apart is not guaranteed to increase the likelihood of polarisation, as shown by the results with the BEBA model.Although the Backfire Effect does move opinions in opposite directions, polarisation is a very unlikely outcome.On the contrary, compared to the FJ and CB models, the BEBA model predominantly produces perfect consensus.In fact, polarisation not only requires opinions to be significantly distant, but also requires the formation of two distant groups including a very similar number of agents, which rarely happens when the initial opinions are random.If the initial opinions have a significant bias, then moving opinions to opposite extremes simply results in a greater majority of agents having one extreme opinion, leading to perfect consensus.
Like the FJ model, the CB model has a significant potential of producing opinion distributions that do not form perfect consensus.Our results for this model suggest that it has a particularly complex behaviour, since the QOL Tables do not display many clearly distinguishable regions associated with qualitative behaviours.However, the average radical trait is shown to have a significant effect on the model outcome.
The analysis of the BC model showed that, although this model is widely related to clustering behaviour, this outcome is not the most likely, and in fact occurs in very particular circumstances, when the initial opinions are polarised or evenly distributed, and the distribution of agent's confidence radii has mean around 0.5 and variance around 0.8.In general, for this model, a uniformly distributed agent parameter set creates more variety in the predicted opinions, because heterogeneous confidence radii result in a network that is not necessarily strongly connected and therefore allows agents to have different opinions.On the other hand, if all the confidence radii are similar, then the agents are either all isolated or strongly connected.
The Agent Parameter ordering plots in Fig. 3 suggest that the mean of the agent parameters is the most significant global parameter in the qualitative description of the model behaviour.On the one hand, this is not surprising, as the mean of a set of numbers is one of its most fundamental properties.In this context, the mean parameter characterises the average agent and is thus a good measure of how the population behaves collectively.For instance, for the FJ model, knowing that the average agent has a high susceptibility is more relevant than knowing that all agents have the same-unknown-susceptibility.On the other hand, it is remarkable that a simple indicator like the mean can capture the overall behaviour of such complex systems.In fact, this is not always the case: there could be models where the Agent Parameter ordering can vary significantly along both axes.The BEBA and BC plots in Fig. 3 suggest how these orderings could look like.
It is worth stressing that our analysis does not require any assumption on the convergence to an equilibrium.Our proposed QOL Tables and QOL Figures provide a static snapshot of the evolution of social systems, suitably sampling the parameter space and the space of possible initial conditions, after a fixed number of time steps.If, for the considered model and the chosen simulation horizon, an equilibrium has been reached at the end of the simulation, then QOL Tables and QOL Figures capture the asymptotic behaviour of the system.On the other hand, our probabilistic analysis approach could also be adopted to study non-stationary asymptotic behaviours, which have been observed in some models of opinion dynamics 47,48 .In fact, the analysis of a sequence of QOL Tables and QOL Figures corresponding to different simulation horizons can provide further insight into the overall temporal evolution and properties of a given model, revealing transient or periodic behaviours.
If QOL Tables and QOL Figures are used to analyse the asymptotic model behaviour, all the simulations need to have reached equilibrium.Heuristic approaches can be adopted to properly choose the simulation horizon, so as to guarantee that the results of such an analysis are sound and interpretable.A sequence of QOL Tables and QOL Figures for increasing time horizons can be compared and, once no change is detected, all systems can be assumed to have reached equilibrium.As an alternative, we can adopt a procedure similar to that used to construct QOL Tables and QOL Figures, with the Bernoulli trial being the answer to the question: Has the system reached equilibrium?
The proposed technique is part of a recent trend of simulation-based analysis methodologies and frameworks used to study and compare opinion formation models that, because they include various social and psychological traits, are too complex to be thoroughly assessed analytically.As such, the proposed QOL analysis can be used in conjunction with other recently proposed approaches, such as the Transition Tables 38 , the Agreement Plot 39 , the Global Unifying Frame 41 , and the distributional measures of Bias, Diversity, and Fragmentation 13 .
Importantly, our proposed QOL analysis technique accounts for incomplete information to reflect reality: for social systems, the model parameters and the initial opinions are extremely challenging or even impossible to estimate, and this crucial aspect cannot be overlooked by assuming exact initial conditions and parameter values.
Another advantage is that the systematic application of the proposed technique allows for a comprehensive exploration of the model capabilities from a wide range of initial conditions and model parameters: different initial opinions (not only uniformly random, which is the case in most simulation results), agent parameters (that are frequently assumed to be the same for all agents), and underlying digraphs (not only regular graphs).This, in turn, offers a clearer and broader picture of the opinion distributions that a model can produce.
On the other hand, this thorough analysis requires a significant computational burden: many simulations need to be performed to obtain meaningful results.For instance, creating the QOL Tables in Fig. 4a has required the simulation of approximately (318)(371)(10, 000) = 1.1798 • 10 9 different social systems.Even with code optimi- sation and multi-core computers, this is a considerable effort.Furthermore, the results only provide qualitative information on the system outcome.
In particular, the amount of time T sim needed to produce the QOL Tables depends on the model, the simu- lation horizon, the density of samples in the parameter and initial condition space, and the number of events (number of Bernoulli trials).The product of the simulation horizon K and the single step calculation time T step approximates the time a single model evolution requires, while the product of N e |{O i }||{P j }| is the total number of simulations (with |X| denoting the cardinality of set X). Since all the simulations are independent, the analysis is easily parallelizeable: if N e |{O i }||{P j }| cores are available, then T sim ≈ KT step ; on the other hand, if a single core is available, then T sim ≈ N e |{O i }||{P j }|KT step .In general, if N cores cores are used to compute the different simulations, then T sim ≈ ⌈N e |{O i }||{P j }|/N cores ⌉KT step , where ⌈•⌉ is the ceiling operator.The single step calcula- tion time, T step , depends in turn on the model, the number of agents, the hardware, and the code implementa- tion.However, in practical applications it will be at most of the order of milliseconds (in a 2.3 GHz Intel Core i5 processor, the value of T step is about 0.343 ms, 0.377 ms, 0.182 ms, and 0.015 ms for the BC, BEBA, CB, and FJ models respectively, implemented in Matlab).In a single core case, with the parameters used for the FJ model ( K = 50 , N e = 10, 000 , |{O i }| = 318 , and |{P j }| = 371 ) a total of 5.9 • 10 10 model time-steps are required, which take approximately 8.64 • 10 5 seconds, or about 10 days.The classification algorithm to decide the Bernoulli trial needs to be executed N e |{O i }||{P j }| times, which can add significant time to the overall computations, so that, even with multiple cores available, the simulations can take several days.The dependence of T step on the num- ber of agents varies with the model.For models like the weighted-median model 11 , T step is independent of the number of agents.For models like the BC model or the CB model, the number of operations grows polynomially with the number of agents.In our simulations, we considered 100 agents, a number that is big enough for the system to exhibit complex behaviours, but small enough to make the simulation time reasonable.Although the algorithm is easily parallelizeable, since the simulations are independent, the computational burden is a limiting factor when exploring more dense parameter space and initial condition space, and longer simulation horizons.
The choice of simulation horizon K, number of trials N e , and size of the families {O i } and {P j } depends on the analysis requirements and available resources.The simulation horizon can be chosen by fixing the time interval to be investigated and estimating the number of opinion changes the agents can have in this time interval.For our case, we considered a time interval of 5 years and estimated 10 opinion changes per year (as done in 38 to align with the World Values Survey timeline), resulting in K = 50 .The number of trials N e depends on the desired confidence.As seen in Eq. (1), δ → 0 as N e → ∞ .In our case, we chose N e = 10, 000 so that all the calculated probabilities had an uncertainty of at most 0.01 (see Fig. 1 in the Supplementary Information).Equation (1) can thus be used when determining the needed number of trials.The reasonable size of the families {O i } and {P j } depends on the purpose for which the QOL Tables are constructed.Since we plotted the QOL Tables for visual inspection, the appropriate size of {O i } and {P j } was determined by the desired image resolution.On print, the plots of the QOL Tables would be at about one inch per side, and 300 pixels/inch is a standard high-quality printing resolution: hence, we chose the number of samples so that tables had around 300 rows (size of {O i } ) and columns (size of {P j } ).Given the limited resolution that the human eye can perceive, taking more samples would provide comparable results, while requiring more resources.Of course, the size can be freely chosen when implementing the methodology so as to comply with different requirements.
Another factor that could affect the size of {P j } is the parameter space of the model.The parameter space is bounded in the FJ model (susceptibility lies in the interval [0, 1]), in the BC model (the confidence radius lies in the interval [0, 2]) and in the CB model (the possible agent traits are bounded in the simplex α ≥ 0 , β ≥ 0 , γ ≥ 0 , and α + β + γ = 1 ).Conversely, for the BEBA model the entrenchment can be any positive value ( ρ > 0 ).When the parameter space is bounded, taking a family {P j } that covers all the parameter space is equivalent to carrying out a sensitivity analysis, exploring the range of all possible parameters; therefore, our approach lends itself to be used as a tool for sensitivity analysis, showing the impact of choosing different parameter values.On the other hand, when the parameter space is unbounded, a preliminary sensitivity analysis is needed to determine the portion of the parameter space that {P j } will cover.In the case of the BEBA model, we found out that choosing an entrenchment value higher than 7 yields extreme model behaviours that lead to complete polarisation almost immediately, and hence we considered entrenchment values bounded as 0 < ρ < 7.
The proposed technique is very versatile and can be used to study any agent-based opinion formation model, as long as the opinion of every agent is a bounded real number.Additionally, the Bernoulli trial can be modified to be anything the study requires.So, for instance, instead of assigning the predicted opinion distribution to one of the five categories, the Bernoulli trial could be asking whether the predicted opinion distribution meets certain desired criteria (such as being in a given subinterval).Another possibility is to use the probabilistic analysis to study the effect of the graph topology on the model evolution.Here, we have considered digraphs with a Small-World network topology because they possess high clustering coefficient and low diameter, characteristics that have been observed in real-life societies 49,50 .However, by simply changing the set N to contain digraphs with e.g.Scale-Free or Random topologies, and comparing the resulting QOL Tables and QOL Figures, the overall effect of the network topology on the model behaviour can be studied via our proposed methodology.
The predictions offered by recent agent-based opinion formation models explain some opinion evolutions seen in real life, but they still do not account for other ones.To gain further insight into the complex mechanisms behind opinion formation, increasingly more sophisticated and intricate models are likely to be proposed, and their analysis will require thorough numerical techniques such as the one we proposed in this paper.

( 4 )
P B = P B O, P = P B x, |x|, , σ ( ) (a) The location of Oi in the Agreement Plot shows the average and the average of the absolute values of the elements in the set.(b) The location of Pj in the Parameter Plot represents the average and the variance of the susceptibility of the elements in the set.

Figure 1 .
Figure 1.Sets O i and P j used in the QOL analysis of the Friedkin-Johnsen model.Taking a point in panel (a), with coordinates (x, |x|) , and a point in panel (b), with coordinates ( , σ ( )) , we can assess the probability P B that a population starting from an initial opinion distribution with average x and average of the absolute values |x| , whose susceptibility has average and variance σ ( ) , produces a predicted opinion distribution categorised as B , with B ∈ {PC, Co, Po, Cl, Di}.

Figure 2 .
Figure 2. Construction of QOL Tables.Each of the points in the Agreement Plot has coordinates (x) i , (|x|) i corresponding to a set O i : they are ordered as O 1 , O 2 , O 3 following the red arrows.Analogously, each of the points in the Parameter Plot has coordinates (p a,j , p b,j) corresponding to a set P j : they are ordered as P 1 , P 2 , P 3 , P 4 following the red arrows.Then, the QOL Tables are constructed, whose rows are associated with points in the Agreement Plot and whose columns are associated with points in the Parameter Plot, in the corresponding order.Therefore, the cell in position (i, j) in the QOL Tables is associated with point O i = (x) i , (|x|) i in the Agreement Plot and with point P j = (p a,j , p b,j ) in the Parameter Plot, and contains a value between 0 and 1 that represents the probability P B .
(d) BC model, with agent parameter sets r (confidence radius) characterised by average r and variance σ(r).

Figure 3 .
Figure 3. Ordering of the points in the Agreement Plot and in the Parameter Plot for the four considered models.Each point (x, |x|) in the Agreement Plot corresponds to a set O of initial opinion distributions.Each point (p a , p b) in the Parameter Plot is related to a set P of agent parameters.The order of the points, indicated by the colorbars, is used when plotting the corresponding QOL Tables as explained in Fig.2: Agreement Plot points are associated with rows, Parameter Plot points with columns of the QOL Tables.This leads to an intuitive representation: for instance, cells near the bottom left corner of the QOL Tables correspond to high rows and low columns, which for the FJ model, according to the order in Fig.3a, are associated with societies with high x and low |x| (i.e., initially polarised), and with high and low σ ( ) (i.e., very susceptible).

Figure 4 .
Figure 4. QOL Tables and corresponding histograms for the FJ, CB, BEBA, and BC models.First row: QOL Tables for B ∈ {PC, Co, Po, Cl, Di} .The rows corresponds to points (x, |x|) in the Agreement Plot and the columns to points (p a , p b ) in the Parameter Plot, both ordered as in Fig.3: each cell is associated with a set of initial opinion distributions and a set of agent parameter sets.The cell colour represents the probability P B , computed as per Eq. 1 (with 10000 samples), that the predicted opinion distribution after 50 time steps is categorised as B ∈ {PC, Co, Po, Cl, Di} .For details on the construction of QOL Tables, see Fig.2.Second row: histograms corresponding to the QOL Tables, visualising the number of cells in the table whose probability value P B is in a given interval; the bar colour reflects the corresponding probability.

Figure 5
Figure 5 presents the QOL Figures for the FJ, CB, BEBA, and BC models.

Table 1 .
Probability of obtaining an opinion distribution of the corresponding qualitative category (Perfect consensus, Consensus, Polarisation, Clustering, Dissensus) for each opinion formation model (Friedkin-Johnsen, Classification-based, Backfire Effect and Biased Assimilation, and Bounded Confidence) after 50 time steps.The cell colour saturation is proportional to the reported probability.The number of samples used to compute this probability, for each model is: 1198860000, 1154340000, 1014420000, and 1179780000 respectively (10,000 times the number of cells in the corresponding QOL Tables).Add all the Bernoulli trial successes of simulations for which the initial opinion distribution belonged to O i .Numerically, this is equivalent to: Vol:.(1234567890) Scientific Reports | (2023) 13:20152 | https://doi.org/10.1038/s41598-023-46789-3www.nature.com/scientificreports/(i)